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Unconventional symmetries of the order parameter allowed some researchers to maintain that 
a purely repulsive interaction between electrons provides superconductivity without phonons in a 
number of high-temperature superconductors. It is shown that the Cooper pairing in p and d states 
is not possible with the realistic Coulomb repulsion between fermions at relevant temperatures in 
any dimension. 
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In the theoretical analysis, the pairing mechanism of 
carriers could be not only phononic as in the BCS the- 
ory or its strong-coupling bipolaronic extension but 
also excitonic, plasmonic, magnetic, kinetic, or due to 
some purely repulsive interaction combined with the un- 
conventional pairing symmetry of the order parameter. 
Actually, following the original proposal by P. W. An- 
derson, many authors Q assumed that the electron- 
electron interaction in high-temperature superconduc- 
tors was strong but repulsive providing high T c without 
phonons via superexchange and/or spin- fluctuations in 
the d-wave pairing channel (1 = 2). A motivation for this 
concept can be found in the earlier work by Kohn and 
Luttinger (KL) Q, who showed that the Cooper pairing 
of fermions with any weak repulsion was possible since 
the two-particle interaction induced by many-body ef- 
fects is attractive for pairs with large orbital momenta, 
Z> 1. 

While the KL work did not provide the specification 
of the actual angular momentum of condensed Cooper 
pairs, Fay and Layzer [|| found that a system of hard- 
sphere fermions condenses at low densities into a p- 
orbital state (I — 1). The critical transition tempera- 
ture T c of repulsive fermions was estimated well below 
0.1 K with very little enhancement due to a flatness of 
the Fermi surface [5j- In two dimensions (2D) the KL ef- 
fect is absent for the parabolic band dispersion, but the 
d-wave pairing appears with the repulsive Hubbard U 
potential when tight-binding corrections to the electron 
energy spectrum are taken into account @ . A few subse- 
quent studies argued that the KL pairing with moderate 
values of angular momenta (p or d) was impossible for 
charged fermions with the realistic finite-range Coulomb 
repulsion in the dense limit, r s < 1, where the perturba- 
tion expansion is justified p-0|- Nevertheless, formally 
extending the perturbation results to the opposite dilute 
limit, r s — > oo, Ref. @ found the p-wave pairing with T c 
about 0.01 K. 

More recent studies claimed that weak repulsive in- 
teractions combined with lattice-induced band-structure 



effects do result in higher values of T c in a spin singlet 
d-wave channel near half- filling " encouragingly similar to 
what is found in the cuprate high-temperature supercon- 
ductors" Q. More surprisingly these studies claim that 
"in the small r s limit, the Coulomb interactions are suf- 
ficiently well screened that it may be reasonable to treat 
them as weak and short-ranged" @. 

Here we show that the weak and intermediate Coulomb 
repulsions do not result in the Cooper instability in 
p- or d-wave channels in 3D and 2D with or without 
lattice-induced band-structure effects. Based on some 
variational simulations we also argue that the repulsive 
hard-sphere (i.e. Hubbard U) model does not account 
for high-temperature superconductivity in the strong- 
coupling regime either. 

In the framework of the BCS theory the symmetry of 
the order parameter A(p) is found from the linearised 
integral equation 

with the most negative eigenvalue A. Here the integral 
is taken over the Fermi surface (or contour in 2D), is 
the normalisation volume in 3D (d = 3) or area in 2D 
(d = 2), vf is the Fermi speed, and K(p,p') is the two- 
particle vertex. The sum of the first and second order 
diagrams. Fig. (HJ), for the two-particle vertex is evaluated 
from HQ: 

K( P ,p') = Kp-pO + ^p-pO^PMp-pO 

k 

- V (k + p')-«(k-p)]Q(p-p',k) 
- 5>(p-kMk + p')Q(p + p',k) (2) 
k 

with Q(q,k) = (/ k - / k _ q )/(£ k - ^k-q)' 4 is the 
Fermi-Dirac distribution function, and E-^ is the Bloch- 
band dispersion. 

The second-order terms in Eq. ([2]) are prohibitively dif- 
ficult to evaluate when the Fourier transform of the re- 
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pulsive potential, w(q) depends on q, so that most previ- 
ous and recent studies [j| Q confined to the hard-sphere 
(Hubbard U) repulsion with u(q)=constant. In this case 
the first order does not contribute to A in any unconven- 
tional channel with I 1, and the only (attractive) con- 
tribution comes from the last sum in Eq.Q providing p- 
wave pairing in 3D for parabolic dispersion [4| and d-wave 
pairing in 2D for the tight-binding dispersion @, H, § . 

Different from those studies we consider here a more 
realistic Coulomb repulsion with u(q) = Aire 2 /Q(q 2 + k 2 ) 
in 3D and v(q) = 2ire 2 /£l(q + k) in 2D, where k is the 
inverse screening length. To elucidate the role of a finite 
potential radius we first take k as an independent large 
parameter k» kp, where Hkp is the Fermi momentum. 
Such a short screening length could be due to another 
component of heavy carriers. In this limit, which mimics 
the Hubbard U model, one can neglect the q dependence 
of the potential in the second-order diagrams, while tak- 
ing it into account in the first-order contribution. Then 
the first sum in Eq.([2]) cancels. In the case of the 3D 
parabolic dispersion = k 2 /2m one can expand the 
order parameter in a series of the Legendrc polynomi- 
als, P;(cos 0) on the spherical Fermi surface [|[ to obtain 
from Eq.([T]) the Z-channel eigenvalue as 
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FIG. 1: (Color online) Angular momentum of the order pa- 
rameter as the function of the dimensionless repulsion s — 
e 2 /irhvF in the 3D Coulomb gas. Inset shows the first order 
(a) and the second-order (b,c,d,e) contributions to the two- 
particle vertex. 
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dQ sin(6)P ; (cos 6)r(cos 6), 
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2(1-cos6) + (k/A:f) 2 

i - cos e 



, V2 + vi + cos e 

In— . (4) 



k 4 4-V2(l + cos6) V2 - VI + cos 9 



where we define the small expansion parameter as s = 
e 2 /(irHvp) pa r s /6 (r s is the dimensionless Wigner-Seitz 
radius, and vp = hkp/m is the Fermi speed). 

Integrating in Eq.Q over the scattering angle yields 
A x /s = (fc F / K ) 4 [4/3 - 2s(21n2 - l)/5] and A 2 /s = 
{k F /K) 4 [16(kp/n) 2 /15 - 4s(8 - llln2)/105] for p and 
d -channels, respectively, at strong screening, njkp 3> 1. 
For the p-wave symmetry a negative A appears only if 
s > 10/[3(21n2 — 1)] pa 9 no matter how small the 
screening length is. This surprising result is due to a 
nonvanishing repulsive contribution to the p-wave chan- 
nel of the lowest-order q-correction to the Hubbard repul- 
sion (oc q 2 / k 4 ), which is proportional to the same power 
in the screening length (k~ 4 ) as the attraction. Dif- 
ferent from the p-channel this lowest-order q-correction 
does not contribute to the repulsion in the d-channel, 
so that A2 becomes negative at some critical screening, 
{n/kp) 2 > 420/[15s(8- 11 In 2)] pa 75/s. However, it is 
unrealistic to find the screening length one order of mag- 
nitude smaller than the lattice constant in real solids, so 
that neither p- nor d-wave pairing can be realised in the 
3D Coulomb gas in the weak- and intermediate-coupling 
regimes, where s < 1, if the screening length, ft" 1 , is 
chosen to be small. 



Actually the screening length in the single-component 
dense Coulomb gas, where the perturbation expansion in 
powers of s makes sense, is large, (n/kp) 2 = 4s < 1, so 
that it is unreasonable to treat the Coulomb repulsion as 
weak and at the same time as short-ranged, contrary to 
the suggestion of Ref.@. Following Refs.fl, 0] we now 
include all buble diagrams in the screened potential re- 
placing the screening momentum k by the static Lind- 
hard function, 
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k 2 F ~q 2 /4 ^kp + q/2 



2qk t 



k F - q/2 



(5) 

It is sufficient to use the Fourier component of the 
potential with the long-wave screening, s(0) = 4s and 
the backward scattering p' = — p with p = q_/2 where 
q_ = p p' in diagrams of Fig.(p]c,d), and the forward 
scattering with q + = p + p' pa 2p in the diagram Fig.([T] 
e) Q. Then the all diagrams are readily simplified with 
the following result for T(cos 8) = T a + T b + T c + T d + T e 
in the eigenvalue equation (J3J 



r c 

In 



r fc = 



*(<?-)' 

4s 2 



kdk 



[q 2 -+s(q-)]q- Jo (k 2 - q 2 _ / A + As) 
q-/2 + k k 2 + q 2 _/A + As + kq^ 
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FIG. 2: (Color online) The absolute value of the coupling con- 
stant for the 2D Hubbard model with the nearest-neighbor 
hopping as the function of the chemical potential for four un- 
conventional symmetries of the order parameter. Inset shows 
the angular dependence of the d x 2_ y 2 state similar to that 
found in Ref. M. 



FIG. 3: (Color online) U — fi phase diagram of 2D electrons 
on the square lattice showing no d- or p-wave symmetry of 
the ground state for weak repulsion U. There is a critical 
value of U at any doping below which the ground state is not 
superconducting (NSC). 



where q 2 . = 2(1 ± cos 6). 

Eqs. (|6l3p allows us to find the symmetry (i.e. the an- 
gular momentum I) of the order parameter with the most 
negative value of A, shown in Fig.([T]) as the function of 
the interaction strength. As in the case of independent 
screening we do not observe the p and d-wave pairing 
in the whole region of the Kohn-Luttinger perturbation 
expansion and beyond (up to s = 3) in agreement with 
Refs. [1, 0] • There is a pairing in higher momentum states, 
1^3, but the corresponding eigenvalues are numerically 
so small (A3 « 0.0011 for s = 3), that the correspond- 
ing T c « {EpfUs) exp(— 1/A) is virtually zero for any 
realistic Fermi energy Ep. 

Finally, let us analyse the unconventional pairing in a 
two-dimensional electron gas on the square lattice with 
a tight-binding energy dispersion, 

Ep = —2t[cos(p x a/h) + cos(p y a/h)} — [i (7) 

(a is the lattice constant). For the half- filled band the 
Fermi level, /i is found at the van-Hove singularity (vHs) 
of the density of states, \x — 0, so that one might expect a 
strong enhancement of the unconventional T c near half- 
filling due to vHs proximity Q. The repulsion between 
electrons is modeled by a strongly screened Coulomb po- 
tential with the Fourier component 27re 2 /£7(<7 + k) « 
27re 2 (l — q/n)/QK, where the inverse screening length is 
taken as an independent large parameter, na 1. 

Similar to the 3D case, the short screening length al- 
lows one to neglect q dependence of the potential in 
the second-order diagrams slightly overestimating their 
contribution, while taking q into account in the first- 
order contribution. Then the diagrams b,c,d in Fig.([T]) 



cancel each other. The remaining second-order contri- 
bution (e), proportional to the static susceptibility of 
the noninteracting electrons, can be reduced to a one- 
fold integration, so that the two-particle vertex con- 
tribution to the unconventional pairing is expressed as 
K[p,p') = (2Tre 2 /(flaK 2 ))T(ap/h,ap'/h) with 



r(k,k') = -|k-k'|- 
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( ai +b)(a 2 -b) 
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Ol,2 = 

sin^^). 
z = cos 



sin(a: + q y /2) sin(q y /2) tan(±z/2 + q x /A) + 
b = [sin(^/2) 2 - sin(a; + q y /2) 2 sm(q y /2) 2 ] 1 / 2 , 
^ 1 [-cos(a;) - /2)], p, = n/2t, q = k + k', 



and k m = cos (—1 — //). Here one assumes that 
I sin(g 2; /2) | > | sin(<7 y /2) | and if otherwise, one should 
replace q x ^ q y . No matter what the screening length is, 
the relative contribution of the second-order diagram (e) 
depends on a single dimensionless interaction parameter 
U = e 2 /at r s , which is supposed to be small in the 
framework of the KL approach. 

Using Eq.® we can solve the eigenvalue problem, 
Eq.([T]), numerically as in Ref. Q by discretization of 
the Fermi surface and diagonalization of the kernel in 
the following integral equation, 



AA(k) 



U 



47r(fta) 



dk 



,r(k,k')A(k') + r(k,-k')A(-k') 



jsin(fc')sin(fc !y )|i/2 



(9) 

Here vectors k and k' are taken on the Fermi surface, so 
that k' y in the integral is defined via k' x using cos(7rfc( / ) + 
cos(7rfc^) = — fi. 

With 800 discretization points we reproduce fairly well 
the results of Ref. 1 81 for the 2D Hubbard model, if we 
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FIG. 4: (Color online) The absolute value of the coupling as 
the function of the chemical potential for two unconventional 
order parameters. 



drop the first term in r(k, k'), Eq.© arising from the 
expansion of the bare Coulomb potential in powers of 
1/k, Fig.©. As shown in the inset of Fig.© the ground 
state is B\ g spin singlet with the d-wave symmetry x 2 —y 2 
close to the half-filling for any weak Hubbard repulsion. 

In fact, there is no reason to neglect the first-order 
q-correction in the two-particle vertex Eq.©, since this 
correction is larger than the spin-fluctuation contribution 
at any screening, if the repulsion is weak. Using the cor- 
rect vertex Eq.© qualitatively changes the ground state, 
Fig.©. Contrary to Ref. [1,1111] neither p- nor d-wave 
pairing are possible in the ground state at any filling until 
the effective interaction becomes so strong that the per- 
turbation theory does not apply, U > 1, Fig.©, similar 
to the 3D case discussed above. As in the continuum 3D 
Coulomb gas, Fig.©, there is a pairing in higher momen- 
tum states (e.g. A2 g with the symmetry (x 2 — y 2 )xy), but 
as shown in Fig.© the corresponding eigenvalues are nu- 



merically very small, and the corresponding T c is about 
zero. 

In cuprate superconductors and many other metal- 
lic compounds the Coulomb interaction is rather strong 
r s 3> 1, so that the perturbative KL approach might have 
no direct relevance to these materials. Different numeri- 
cal techniques have been applied to elucidate the ground 
state of the repulsive Hubbard model in the intermedi- 
ate to strong-coupling regime, U > 1 sometimes with 
conflicting conclusions. In particular, recent studies by 
Aimi and Imada [Io| using a sign-problem- free Gaussian- 
basis Monte Carlo (GBMC) algorithm showed that the 
simplest Hubbard model with the nearest-neighbor hop- 
ping has no superconducting condensation energy at op- 
timum doping. This striking result was confirmed in the 
variational Monte Carlo (vMC) studies by Baeriswyl et 
al. [ll|, who found, however, some condensation energy 
away from the optimum doping and also adding next- 
nearest neighbor hoppings. Importantly, a similar vMC 



method [12[ found that even a relatively weak finite-range 
electron-phonon interaction with the BCS coupling con- 
stant A ~ 0.1 induces a d-wave superconducting state 
in strongly correlated metals with the condensation en- 
ergy several times larger than can be obtained with the 
Hubbard repulsion alone. Moreover, the unconventional 
superconductivity has been shown to exist due to a finite- 
range electron-phonon interaction [e| [H| without the 
need for additional mechanisms such as spin fluctuations. 

In conclusion, we have shown that the p- and d-wave 
Cooper pairing from the weak Coulomb repulsion is not 
possible between fermions at any screening length and in 
any dimension. Pairing in higher momentum states (I ^ 
3) has virtually zero T c for any realistic Fermi energy. 
The unconventional pairing from the strong Coulomb re- 
pulsion is not possible either since the corresponding con- 
densation energy, if any, is several times lower than the 
condensation energy caused by the electron-phonon in- 
teraction. 
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